APPROXIMATION OF THE LONG-TERM DYNAMICS 
OF THE DYNAMICAL SYSTEM GENERATED BY 
THE TWO-DIMENSIONAL THERMOHYDRAULICS 

EQUATIONS 

FLORENTINA TONE 

Abstract. Pursuing our work in [TH], ^7\, [30], [S], we consider in 
this article the two-dimensional thermohydraulics equations. We 
discretize these equations in time using the implicit Euler scheme 
and we prove that the global attractors generated by the numerical 
scheme converge to the global attractor of the continuous system 
as the time-step approaches zero. 
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1. Introduction 

In this article we discretize the two-dimensional thermohydraulics 
equations in time using the implicit Euler scheme, and we show that 
global attractors generated by the numerical scheme converge to the 
global attractor of the continuous system as the time-step approaches 
zero. In order to do this, we first prove that the scheme is if^-uniformly 
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stable in time (see Section HI) and then we show that the long-term 
dynamics of the continuous system can be approximated by the discrete 
attractors of the dynamical systems generated by the numerical scheme 
(see Section [5]). 

In the case of the Navier-Stokes equations with Dirichlet boundary 
conditions, the i7^-uniform stability of the fully implicit Euler scheme 
has proven to be rather challenging. However, using techniques based 
on the classical and uniform discrete Gronwall lemmas, we have been 
able to show the if^-stability for all time of the implicit Euler scheme 
for the Navier-Stokes equations with Dirichlet boundary conditions 
(see [20] )• The if ^-stability has also been established. More precisely, 
the if^-stability has first been proven in the simpler case of space pe- 
riodic boundary conditions (see ^7\), and then extended to Dirichlet 
boundary conditions (see [H]); the magnetohydrodynamics equations 
are also considered in [T8] . 

Our first objective in this article is to extend the if^-uniform stability 
proven in [20j for the Navier-Stokes equations with Dirichlet boundary 
conditions, to the thermohydraulics equations. In order to do so, we 
divide our proof into three steps. First, we prove the L^-uniform sta- 
bility of both the discrete velocity f " and the discrete temperature 6'" 
(see Lemma [3^ below) . Then, using techniques based on the classical 
and uniform discrete Gronwall lemmas, we derive the if^-uniform sta- 
bility of f" (see Theorem 14. II below) . which we will use in Subsection 
14.21 in order to establish the if^-uniform stability of 6"" (see Theorem 
14.21 below). Besides the intrinsec interest of considering the thermohy- 
draulics equations, the new technical difficulties which appear here are 
related to the specific treatment of the temperature with the necessary 
utilization of the maximum principle. Furthermore, we have simplified 
some steps of the proof as compared to [20] . 

Our second objective in this article is to employ the technique de- 
veloped in [5] to prove that the global attractors generated by the fully 
implicit Euler scheme converge to the global attractor of the contin- 
uous system as the time-step approaches zero. When discretizing the 
two-dimensional thermohydraulics equations in time using the implicit 
Euler scheme, one can prove the uniqueness of the solution provided 
that the time step is sufficiently small. More precisely, the time re- 
striction depends on the initial value, and thus one cannot define a 
single- valued attractor in the classical sense. This is why we need to 
use the theory of the so-called multi- valued attractors, which we briefiy 
recall in Subsection 15.11 
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2. The thermohydraulics equations 

Let fl = (0, 1) X (0, 1) be the domain occupied by the fluid and let 
62 be the unit upward vertical vector. The thermohydraulics equations 
consist of the coupled system of the equations of fluid and temperature 
in the Boussinesq approximation and they read (see, e.g., |6], [T5]): 

dv 

(2.1) — + ivV)v-uAv + Vp = e2iT-Ti), 

(2.2) ^ + . v)T - kAT = 0, 

(2.3) dwv = 0; 

here v = (f i, f 2) is the velocity, p is the pressure, T is the temperature, 
Ti is the temperature at the top boundary, X2 = I, and z/, k are positive 
constants. We supplement these equations with the initial conditions 

(2.4) v{x,0)=Vo{x), 

(2.5) T(a;,0)=T°(x), 

where fo : f2 — )■ M^, T'^ : — > M are given, and with the boundary 
conditions 

(2.6) V = at X2 = and X2 = 1, 

(2.7) T = Tq = Ti + 1 at X2 = and T = Ti at X2 = 1, 

and 

p, f , T and the first derivatives of v and T are periodic 

f2 8) 

of period 1 in the direction xi, 

meaning that 0|2;i=o = 0|xi=i for the corresponding functions 0. 
Letting 

(2.9) e = T-To + X2, 

and changing p to 



(2.10) p-[x2 2 
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equations f l2.ip -( l23|) together with the boundary conditions f l2.6p -( l2r8l) 
become 

dv 

(2.11) — + {vV)v - u/\v + Vp = 626, 
at 

do 

(2.12) — + {vV)e-V2-KAe = 0, 

(2.13) divw = 0, 

(2.14) V = at X2 = and X2 = 1, 

(2.15) 6 = at X2 = and X2 = I, 

(2.16) ([21]) holds with T replaced by ^. 

These equations are supplemented with the initial conditions 

(2.17) v{x,0) = vo{x), 

(2.18) e{x, 0) = r°(x) -T0 + X2 =: Ooix). 

For the mathematical setting of the problem we define the space H = 
Hi X H2, where 

(2.19) 

Hi = [v e L^{Vtf, div = 0, V2\x2=o = V2\x2=i = 0,t;i|xi=o = ^^iUi=i} , 

(2.20) H2 = L'^in), 

and we denote the scalar products and norms in Hi, H2 and H hj {■,■) 
and I ■ |. 

We also define the space V = Vi x V2, where 
(2.21) 

Vi = {v e H^i^Y, f 1x2=0 = ^^1x2=1 = 0,^1x1=0 = ^1x1=1, div V = 0] , 

(2.22) V2 = {9 E H^{ri), 9\x2=o = 9\x2=i = 0,6'|xi=o = ^Ui=i} • 
The space V2 is a Hilbert space with the scalar product and the norm 

(2.23) ((0, ^)) = /" V0 ■ dx, U\\ = ^M^, 



and we have the Poincare inequality (see, e.g., [T^, page 134) 
(2.24) |0|<||0||, V0eViorl^2. 

We denote both scalar products and norms in Vi and V by ((■,■)) and 

Ml- 

Let D{A) = D{Ai) X D{A2), where 



(2.25) DiA,) = \veV,nH\nf, 



dv 




dv 




dxi 


xi=0 


dxi 


xi=i\ 



1,2, 
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and let A be the hnear operator from D{A) into H and from V into V 
defined by 

(2.26) {Aui, U2) = a{ui, U2), ^ Ui = {vi, Oi] e D{A), i = l,2, 
with 

(2.27) a{ui, U2) = u{{vi, V2)) + ^{{61, 62)). 

We consider the trihnear continuous form b on V, defined by 

(2.28) b{ui,U2,U3) =bi{vi,V2,V3) + 62(^^1, 6^2, 6^3), Vu^ = {fi,6'i} e V, 
where 

(2.29) b,{y,w,z)=Y, m-p;Zjdx,^y,w,zeH\Qy, 

(2.30) b2{y,(jy,ij) = y2 [ Vi^i^ dx^^V ^ H\n)\<f),^ e H\n). 

i=i -^^ * 

The form bi is trihnear continuous on Vi x Vi x Vi and enjoys the 
foUowing properties: 

(2.31) My,w,z)\ < c,|y|V2||^||V2||^|||^|i/2||^||i/2^ Wy,w,ze V„ 

.2 32^ My,w,z)\<c,\ynA,ym\w\\\zl 

^yeD{A,),weV,,zeH,, 

|6i(y,«;,^)|<c,|t/r/^|b||V2||^u;||V2|Ai^|V2|^|^ 
'^yEVuweD{Ai),zEH,, 

(2.34) bi{y,w,w) = 0, ^y,weVi, 
the last equation implying 

(2.35) bi{y,w,z) = -bi{y,z,w), 'iy,w,zeVi. 

The form 62 is trihnear continuous on Vi x V2 x V2 and enjoys the 
following properties, similar to fl2.31l) - fl2.35l) : 

(2.36) \b2{y,^,n < c,|l/^/ly|^/^||0|||v^^/^||^|^/^ vy,^,^ g v2, 

.2 37^ MyAM<c,\ynA2yn\mi 

Vy G DiA2), G F2, ^ e H2, 

i62(?/,0,v^)i<c,iyr/^ibir/^ii0r/v20r/'i^i, 

Vy G ^1, G DiA2), G H2, 
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(2.39) 62(z/,0,0) = O, VyeVi, 6^2, 
the last equation implying 

(2.40) 62(z/,0,V^) = -&2(z/,V^,0), VyG Vi, 0,V^e ^2. 

We associate with b the bilinear continuous operator B from V x V 
into V and from D{A) x -D(A) into H, such that 

(2.41) {B{ui,U2),U3)v',V = b{ui,U2,U3), Wui,U2,U3eV. 

We also define the continuous operator in H 

(2.42) Ru = -{626, V2}, u = {v, e}. 

For more details about the function spaces D{A), V and H, as well as 
the operators A, B, R and b, the reader is referred to, e.g., [T5]. 

In the above notation, the system (I2.1ip - (l2.13p can be written as 
the functional evolution equation 

(2.43) ut + Au + B{u) + Ru = 0, u{0) = uq = {vq, Oq}. 

In the two-dimensional case under consideration, the solution to the 
thermohydraulics equations is known to be smooth for all time (cf. 
[T5]). Using the maximum principle for parabolic equations, one can 
show that e G L°°(M+; L'^i^l)) and the velocity u is bounded uniformly 
for all time by 

(2.44) l^(t)li.(np <e-^>oli.(f,). + ^(l-e-^*)' 

where 6*00 = \0\L'=^{R+;L'2{n))- Furthermore, using techniques based on 
the uniform Gronwall lemma (cf. [I5]), one can bound the solution u 
of fl2:i3|) uniformly in V for alH > 0. 

In this article we discretize (I2.43P in time using the fully implicit 
Euler scheme, and define recursively the elements = {v^, of V 
as follows: 

M° = {v^,9^}, where v^(x) = vo(x), and 

(2.45) 

9 (x) = 9o{x) := T (x) — Tq + X2 are given; 

then when u° = {v^, 9^}, ■ ■ ■ , u^~^ = {f 9""^^} are known, we define 
= {v", 9""} eV such that 

(2.46) 



+ ^n,^;) = (62^,^;), G Vi, 

k 



(2.47) 



hen _ QU-l^ ^ ^^^QU^ ^ ^^^^n^ gn^ _ ^ = Q, V^^ G ^2- 

k 
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The above system is very similar to the stationary Navier-Stokes 
equations and the existence of solutions is proven e.g. by the Galerkin 
method, as in [16]. Uniqueness can also be derived as in [16] under 
some conditions. Let us explain this point, which somehow motivates 
the developments in Section [51 For that, we rewrite the system fl2.45p - 
(IZiTj) in the form 

(2.48) 

(t;", v) + i/A;((t;", v)) + A;6i(i;", w", v) - A;(e2^", v) = (t;"-\ v),^v e Vi, 
(2.49) 

(r, 6) + kA;((^", 6)) + fc62(^^", e) - e) = (^"-\ 9),y9 e V2, 

and assume that {w", 6'"} and {■u", are two solutions corresponding 
to the same initial data {vo,9q} G V. Setting v"' = v"' — and 
Qn _ Qn _ Qn^ obtain that {■u", 6'"} is a solution to the following 
system: 

(2.50) 

(5",t;) + z/A;((5",t;)) + A;6i({}",t;",t;) + fc6i(w",'D",w) - k{e2t,v) = 0, Vt; G l^i, 
(2.51) 

{r, 9) + kA;((^", 0)) + A;62(^", ^) + kb2{v'', 9) - k{v^, 9) = 0. V0 G V2, 

Taking t> = in (12.501) and using (I2.34p . we obtain 

(2.52) + i/A;||5"f + A:6i(5", w", 5") - A;(e2^", ^") = 0. 

Using property ( 12.311) of the trilinear form hi and the bound (14.661) 
below on we obtain (for k < kt{\\{vq,9q}\\), with /t7(||{fo, 6'o}||) 

given in Theorem 14.31 below) : 

kbi{v'',v'',v'') < CfeA;|{}"|||{}"||||t;"|| < c^i^sA^I^"! 
^'•''^ <'ik\\vr+'^Kik\v-\^. 
We also have 

A;(e2^",^") < A;|e2^"||i)"| < 

(2.54) I, 1 . 

- 4 II II ;y I I 

Relations (I532D-(EMD imply 

(2.55) ( 1 - -ir.^A;') + -A;||5"f < i/tirP. 
\ u J 2 z/ 

Now taking 9 = 9'"' in (I2.5ip and using (I2.39p . we obtain 

(2.56) |r |2 + Kk\\rf + fc62(^", n - ^^(^^2 , n = 0. 
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Using property (12.361) of the trihnear form 62 and the bound (14.661) 
below on ||6'"'||, we obtain 

(2.57) 

— ^ II II ^ II II oil oil 

We also have 

(2.58) k{v^,9'') < k\v^\\9''\ < A;|{)"|||r|| < -k\\ff' 



4 



Relations (1^3B]) -( E3g]) yield 



(2.59) (1 - cKlk)\r\^ + '^k\\rf < ^A;||{}"ii2 
Adding relations (I2.55P and (12. 59 p . we obtain 

(2.60) 













-k\ij^\^ 




K 



1 _ '^Klk - cKlk - hA + (1 - cKlk - ^k^ '^"'2 



4 2 
Assuming k is sufficiently small, that is 

(2.61) A;<minL7(||H,^o}||), 



2 {^Kl + cKl + i) ' 2 {cKl + I) 

relation (I^IBUD implies = 6^" = 0. Hence, the system (ITia - (ITi7D 
possesses a unique solution, provided that the time-step satisfies the 
constraint (I2.6ip . This is enough to uniquely define the sequence 
{v"', ^"'} for k small enough, but the dependence of the time step k 
on the initial data prevents us from defining a single-valued attrac- 
tor in the classical sense, and this is why we need the theory of the 
multi- valued attractors, that we discuss in Subsection 5.1. 

Our next aims are to prove that the solution = {f 6*"} to the 
discrete system fl2.45p - fl2.47p is uniformly bounded in the \^-norm and 
then to show that the global attractors generated by the numerical 
scheme (12.450 - 02.470 converge to the global attractor of the continuous 
system as the time-step approaches zero. 

In this article we only consider time discretization, we do not con- 
sider space discretization. Important background information on space 
discretization and on various computational methods can be found in 
some of the books and articles available in the literature. On finite el- 
ements, see, e.g., [7j, [S]; on finite differences and finite elements, [TU], 
on spectral methods, [3], [S]- 
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3. //-Uniform Boundedness of t;" and 6^" 

In proving the if-uniform boundedness of f " and 6'", we need first to 
prove a variant of tlie maximum principle for 9^. In order to do so, we 
introduce the following truncation operators (cf. fL5]), that associate 
with the function the functions and given by 

(3.1) V+i.^) = max(</)(x), 0), V-{x) = max(— y9(x), 0). 

Note that, with this notation, we have ip = (f^ — if-, the absolute value 
\ip\ of (f is ip^ + ip- and f+f- = 0. Using these operators, we can prove 
the following preliminary lemma 

Lemma 3.1. If'p,ip ^ L'^i^), then 

(3.2) 2{ip - </.+) > - + - V'+r, 

(3.3) - 2(<^ - ^, ip_) > - + - ^_|2. 
Proof. We have 



(3.4) 



since ip-'p+ > 0. The proof is similar for (13. 3p and the lemma is 
proved. □ 

We are now able to prove the following variant of the maximum 
principle for 6"". 

Lemma 3.2. If v"" and satisfy flTi^ and (I237D; ihen 

(3.5) r = r + r, 

with 

(3.6) X2 - 1 < ^" < X2, 

(3.7) 1^1 < (|0°| + |0°|)(l + 2fi:A;)-t. 

Moreover, there exists Mi = Mi{\6o\), given in fl3.26p below, such that 

(3.8) |^"| < Mi,Vn > 1. 

Proof. Rewriting ( I2.47P in terms of T" = 6'" + Tq — X2, we find: 
(3.9) 

i(T" - T"-\ T) + /s:((T", T)) + 62(^^", T'^, T) = 0, VT G U2, = 0, n > 1. 





= 2(</^+ - 




= 2(</^+ - 




-2(y._- 
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Replacing T by 2k(T"- — To)+ in the above equation and using (13. 2p . 
we obtain: 

|(T"-To)+p-|(T"-i-To)+r 
+ |(T" - To)+ - (T"-i - To)+p + 2A;k||(T" - To)+f < 0. 

Using the Poincare inequahty f l2.24p . we find 

(3.11) |(T"-To)+p<-|(T"-l-To) + |^ 

a 

where 

(3.12) a = l + 2Kk. 
Using recursively (13.111) . we find 

(3.13) |(T" - To)+|2 < (1 + 2/s:A:)-"|(T° - To)+p. 
Similarly, using (13.31) . we obtain 

(3.14) |(r"-ri)_|2 < (l + 2Kfc)-"|(T°-Ti)_|l 
Setting 

(3.15) J^n^fn^ fn^ ^-^j^ fn ^ ^j.n _ j.^^^ _ ^j.n _ 

we find that f"^ = T'^ - (T" - ro)+ + (T" - Ti)_, so that T" = Ti, for 
T" < Ti, f " = T", for Ti < T" < Tq, and f " = Tq, for T" > Tq; in all 
cases 

(3.16) Ti<T"<ro. 
Rewriting (I3.13p -( l3.15l) in terms of ^, we obtain 

(3.17) |(r - a;2)+P < (1 + 2/s:A:)-"|(e° - 0:2)+^, 

(3.18) |(^" - 0:2 + 1)_P < (1 + 2/s:A;)-"|(^° + 

(3.19) r + To - X2 = + (r - X2)+ - {e^^ -X2 + 1)-. 
Setting 

(3.20) r = (0"-X2)+-(^"-X2 + l)-, 

(3.21) ^" = T"-To + a;2, 
equation (I3.19P becomes 

(3.22) r = r + r. 
By fl3.16p . we have 

(3.23) a;2 - 1 < ^" < X2, 
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and by (Km . (KT7\\ and (KW we derive 

m<\{9--x,U\ + \{9--X2 + lU 
■ <(l + 2/€A;)-t(|eO| + |eO|). 

To complete the proof of the lemma, we note that (13.22^ . ( ]3.23p and 
I^M) yield 

(3.25) |r| < + {\el\ + le^i) (i + 2Kk)-^,yn > i, 

and setting 

(3.26) M,{\9o\) = + \9l\ + \9% 

we obtain conclusion fl3.8p of the lemma. □ 
Corollary 3.1. // 

(3.27) k< —, 

then Bi2{fd,2\VL\^^'^) , the hall in centered at and radius 2\VL\^/'^, is 
an absorbing ball for 6'" in . 

Proof. Indeed, let B be any bounded set in and assume that it is 
included in a ball 5(0, R) of L^. It is easy to deduce from (K^ that 
for any 9o G B{0,R), 

(3.28) |^"| < + 2R{1 + 2Kk)-^,^n > 1, 

and using assumption fl3.27p on k and the fact that 1+x > exp(x/2) if a; G 



2 In 



(0,1), we obtain that there exists NQ{R,k) := — — - such that 
61" G 5^2(0, 2|(]|i/2),Vn > N^. This completes the proof of the corol- 
lary. □ 

We are now able to prove the if-uniform boundedness of f". More 
precisely, we have the following: 

Lemma 3.3. Let {f 9"'} be the solution of the numerical scheme 
(12. 46p - (12.471) . Then for every k > 0, we have 

(3.29) < (1 + z/A;)-" \vo\^ + h _ (1 + uk)--] , Vn > 0. 

Moreover, there exists Ki = Ki{\vo\, \9o\), such that 

(3.30) < Ki, Vn > 0, 
and 



m ^ m 

(3.31) vk^Wv^W" < \v''^\^ + -kJ2\9'\' 



Vi = 1, ■ ■ ■ , m, 



V 

j=i j=i 
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m ^ m 

(3.32) Kfc^ll^i^ < |^.-i|2^_;.^|^,|2^ Vi = l, 



Proof. Taking v to be 2kv"' in fl2.46p and using the relation 

(3.33) 2(<^-^,<^) = |<^|2-|^|2 + |<^-^|2, 

as well as the skew property 02.341) . we obtain 

(3.34) - + - t;""^!^ + 2vk = 2A;(e2^'', w''). 

Using the Cauchy-Schwarz inequality and the Poincare inequality (I2.24p . 
we majorize the right-hand side of (13.341) by 



2A;(e2r,w") < 2fc|e2^"| < 2k\9''\\v' 

< 2A;|r|||t;"|| < vk\\v''f + -k\d'^ 



(^■^^) ^ ^, i/inn, n„ . , „ .„||2 , ^ r |fln|2 



V 

Relations (I3.34p and (I3.35P imply 

(3.36) - + - i;"-^^ + z/A; ||t;"f < ^A; l^^p. 
Using again the Poincare inequality (I2.24p . we find 

(3.37) W"? <-\v^-\^ \ —kW? , 

a au 

where 

(3.38) a = l + uk. 
Using recursively (13.371) . we find 



1 1 " 1 



j=l 



(3.39) 

< (1 + z/fc)-^ |t;T + ^ [1 - (1 + ^fc)-"l 



which proves (I3.29p . 

Taking K\ = + ^ relation (I3.30p follows right away. 

Adding inequalities (I3.36P with n from z to m we obtain (13.311) . 

Now, replacing 6 by 2k6"' in (I2.47P and using the skew property 
flOOj) . we obtain 

(3.40) \9'f - + - + = 2k{v^, r). 



LONG-TERM DYNAMICS OF 2D THERMOHYDRAULICS EQUATIONS 13 

Using again the Cauchy-Schwarz inequality and the Poincare inequahty 
fl2.24p . we majorize the right-hand side of f l3.40p by 



(3.41) 



2A;K,r) < 2k\vl^\\9''\ < 2k\v^\\\9'- 



K 



Relations ( KM and flCTl) imply 



(3.42) 



\n\2 



-1|2 



+ Kk\\e''\Y < -k\v 

K 



.n\2 



Summing inequalities (13.42^ with n from i to m we obtain ( 13.32p . □ 
Corollary 3.2. Let 

(3.43) 



I 1 1 

k < mm i — , — 

2k V 



and set po = 2|fi|^/^ + ^^^^ ^ . Then Bh{0,Po), the ball in H centered 
at and radius po, is an absorbing ball for {f 6'"} in H. 

Proof. Let B be any bounded set in H and assume that it is included 
in a ball B{0, R) of H. For any initial data G B, Corollary O 

implies that 

(3.44) < 2|fi|^/^Vn > N^{R,k), 

and then fl3.37p becomes 

1 



< -\v 
a 



n-l\2 



+ 



av 



Vn > Nl{R, k), 



(3.45) 
where 

(3.46) a = l + uk. 

Iterating the above inequality, we find (for any n > Nq{R, k)) 



..n\2 



(3.47) 



1 4 """^^ 1 

i=l 

= (1 + z/A;)-("-^«) |t;^o |2 + [l - (1 + uk)- 



(n-iVl) 



< (1 + uk) 



-(n-7Vl) 



R^ + —{\Q\+2R' 



4 , , 



(by (I329D and KJM) . 
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and using assumption f l3.43p on k and the fact that 1+x > exp(x/2) if x E 
(0, 1), we obtain that there exists Nq{R, k), 

(3.48) N^{R,k) := — \n ^ ^^^u |-r 



uk 

such that < y5|l]|i/Vi/, Vn > A/'qI + A^-^ =: No{R, k). 

We, therefore, have that {w'^,^} G 5^(0, po), for all n > No{R,k), 
which completes the proof of the corollary. □ 



4. V-Uniform Boundedness of v"- and 6"- 

We now seek to obtain uniform bounds for f " and 6^ in V, similar 
to those we have already obtained in H (see fl3.30p and (13. 8 p above). 
In order to do this, we first derive bounds for f " valid on any finite 
interval of time (see Proposition 14. ll below). and then we repeatedly use 
them on successive intervals of time together with (a discrete uniform 
Gronwall) LemmaS3]to arrive at the desired uniform bounds. Once we 
have obtained the ^-uniform bounds on u we can use those, together 
with a new version of the discrete uniform Gronwall lemma, to derive 
the l^- uniform boundedness of 6'". 

4.1. //^-Uniform Boundedness of t>". 
Lemma 4.1. For every k > 0, we have 



(4.1) ||t;"f < K2\\v''-^\\^ + 4^1, Vn > 1, 



where K2 = 2(1 + 2clKf/u^). 

Proof. Replacing v by 2k{v'"' — v"^'^) in f l2.46p . we obtain 

2|t;" - + z/A;||t;"f - yk\\v''-^f + z/A;||t;" - t;"-^^ 

+ 2k 61 (t;", f t;" - t;"-^) = 2k (ea^", t;" - v""-^] 



(4.2) 



Using properties fl2.34p . (12.351) and (I2.3ip of the trilinear form hi and 
recalling (13. 30 p . we bound the nonlinear term as 



(4.3) 



2fc6i(w",^;",w"-w"-i) =2fc6i(t;",t;"-\t;") (by flCTj) . flOSD ) 

< 2c^k\v''\\\v''\\\\v''-^\\ (by flCTD ) 

2 Z-^ 
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We bound the right-hand side of (14. 2 p using Cauchy-Schwarz' inequal- 
ity, fl^:^ and ( KE\\ : 



(4.4) 



- 2 " " z/ ^ 



Gathering relations (14. 2 p through (14. 4p . we find 

„ z^. „ „„o / 2c? 

'2\v" - v" 'I" + 

(4.5) 



2 " " - z/ ^' 

We thus obtain 

(4.6) < ^2||^"~1^ + 

which is exactly conclusion (14. ip of the lemma. □ 
Lemma 4.2. For every k > 0, we have 

(4.7) ciKlk\\v''\\^ - \\v''f + + ^kM^ > 0, Vn > 1, 
where Ci = 27c^/(2z/3). 

Proof. Replacing v by 2A;Aif" in (12.460 . we obtain 

^^'^^ +2z/A;|Ait;"|2 = 2A:(e2^",Ait;"). 



Using property (I2.32p of the trilinear form bi and recalling (I3.30p . we 
have the following bound of the nonlinear term, 

2A;6i(t;",t;",Ait;") <2cbk\v''\'/^\\v''\\\A,v''\'^/^ 

(4-9) J, 27c4 

<^k\A,v^\^ + ^Kfk\\v^\\\ 

Using the Cauchy-Schwarz inequality and recalling (13. 8p . we bound the 
right-hand side of (14. 8 p by 



2A;(e2r,Ait;") < 2A;|r | 

< -klAiv^'l'^ + -kM't. 



16 TONE 

Relations (l4:8|) -( g:T0l) imply 

||w"f - Wv'^-^f + lit;" - w"-if + z/A;|Aiw"|2 

(4-11) OTr-^ 2 

<^Klk\\vX + -kMl 

from which we obtain conclusion (14 of Lemma 14.21 □ 

In what follows, we will make use of the following discrete Gronwall 
lemma, whose proof can be found in |13], [20] . 



Lemma 4.3. Given k > 0, an integer n^, > and positive sequences 
^n, Vn o-nd (n such that 

(4.12) Cn < ^n-iil + kr]n~i) + k(n, for u = 1, ■ ■ ■ , u^, 

we have, for any n G {2, ■ ■ ■ , ra*}, 

^ 5Z '^M + ^ 5Z Ci exp ( A; ^ j + fcCn • 

i=0 ^ i=\ ^ j=i ^ 

Proposition 4.1 (Estimates on a finite interval). Let T > be fixed 
and let K^{-, - , ■) he the function, monotonically increasing in all its 
arguments, given in (I4.29P below. If the timestep k is such that 

(4.14) k < min{/€i, /^sd^^ol, l^ol), ^siWv'l \0o\,T)}, 
where Ki is given by (I3.43p . 

(4.15) «:2(|t^o|,|^o|) 



Z/2 



mciKlMl ' 
(4.16) K,{\\vWe,lT) ^ 



lQc^KlK2Kl{\\v%%lTy 
then 

(u) 

,n||2 _^ _^ 

^ ^ ,2- 



CiKik[K2\\v''f + +^Mi] < -, Vn = 0,--- ,A^:= [T/k\., 



(Hi) 

14 



k"ll'< lk"-'ll' 



l + 2ciKlk (^\v''-^f+'^kMl 
= I,-- - := [T/k\. 



5z/ 
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Proof. We will use induction on n. The induction consists in showing 
that (i) holds at n = 0, and then showing that if (i) holds for n < m — 1, 
then (i) holds for n = m. If n = 0, one can easily see using the definition 
f l4.29p of that (i) is true. Now assume that (i) holds for n < m — 1. 
Then, by f l4.14p . we have that (ii) is true for n < m — 1 and 



(4.17) A.,, = 1 - AciK^k 



,n||2 



+ -kM( > 0, Vn < m - 1. 



From (14. 7p . we obtain that either 



(4.18) 
or 

(4.19) 



,n+l||2 



> 



2ciKfk 
1 + x/a; 



2ciKfk ■ 

The second alternative is not possible, as it implies, with (14. ip . (i) at 
order n and fl4.14p : 

1 < 1 + v^A^ < 2ciKlk\\v''+'^f 

(4.20) 



Then the first alternative holds and, thus, 



(4.21) 
with 

Since 



2ciKfk\\v 



X = AciKlk 



X 



n+l ||2 



< 1 - vT 



n\\2 



+ -kMt 

V 



X ( X 
; < - 1 + - 

1 + yr^ - 2 V 2 



1 - v^r^ 

we obtain 

2c^Klk\\v''^^f < 2ciKlk M|t;"f 
(4.22) ^ 

• l + 2ciKlk 

which implies, together with (ii) at order n: 



< X < 



+ -kM; 

V 



,n||2 



+ -kMl 

V 



(4.23) \v 



n+l||2 ^ 



1 + 2ciKlk 



n\\2 



+ -kMi 

V 



+ ^^kMl 
5z/ 
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which is exactly (iii) at order n + 1. We rewrite (iii) in the form 

(4.24) ^n<^n-i{l + kr]n-i) + kC, 

with 
(4.25) 



2 \ 14 
vir, r]n = 2ciK't ( Wv^'f + -kMn and ( = —Ml 



V 



5u 



and we have that (14.241) holds for = 1, ■ ■ ■ , m. In order to prove that 
(i) holds for n = m, we use the (discrete Gronwall) Lemma [4.31 and we 
compute the following. If i > 0, using fl3.3ip . we obtain: 



(4.26) 



m—l m—1 / o \ 

kVj = 2ciKlk ( 11^' ir + =^kM^ j 



Kl + -Ml{m-i)k 



Similarly, if z = 0, using again f l3.3ip and recalhng fl4.14p and ( I4.16p . 
we find 



(4.27) 



m— 1 m—l ^ „ 

Y kVj = 2ciKlk Y ( ll^i^ + -kMl 

<- + -ciK^ (k^ + -M^mk 1 . 
5 u \ V ' 



Using fl4.26p we also have 
(4.28) 



Y kC exp ( J2 f^Vj j < -^Ml exp | -CiKl ( KI + -Mlmk 

i=l ^ j=i ' ^ ^ 

By ( 14131) . (I12ZD, and ( lOSl) we obtain 
(4.29) 



mk. 



.m||2 



28 



< llf^lr + —Mfmk exp <^ - + -ciKf Kf + -Mfmk 
\ 5z/ J [5 u \ V 

=:Kl{\vX%\,mk), 

which is exactly (i) for n = m. 

We note that the dependence of K^, ks and K2 on |t;o| and 16*01 is 
through Ki and Mi, but those quantities bound all and |^"|, re- 
spectively. 

This completes the proof of the proposition. □ 
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We now want to extend Proposition 14.11 to obtain uniform bounds 
for ||f"||, for all n > 0. In order to do so, we will repeatedly apply 
Proposition 14.11 on finite intervals of time, considering different initial 
values. Since the bound obtained on each interval is an increas- 
ing function in the corresponding initial value considered (see fl4.29p ). 
bounding uniformly all initial values will provide a uniform bound for 
all n > 0. To do so, we need the following (discrete uniform 

Gronwall) lemma, whose proof can be found in [2D] (see also 



Lemma 4.4. We are given k > 0, positive integers ni,n2,n^, such that 
ni < n^: , ni + n2 + I < n^: , and positive sequences rjn, (n such that 

(4.30) < ^n-i(l + kr]n-i) + k(n, for u = Ui, ■ ■ ■ , n^. 

Assume also that 



n'+n2 n' +n2 

k ^ 'nn< ai{ni,n^), k ^ (n < 02(^1, n=,) 



(4.31) 



n=n' 
n'+n2 



n=n' 



for any n! satisfying ni < n' < — n2. We then have, 
(4.32) a < (^^^^ + «2(ni,n.)) e-(--), 

for any n such that ni + n2 + 1 < n < n^,. 

We are now in a position to give the main result of this section, that 
is, to derive a uniform bound for for all n > 1. 

Theorem 4.1. Let Vq G Vi, 6q E H2, and {v",9'^} be the solution of 
the numerical scheme (12.461) - fl2.47p . Also, let r > 4ki be arbitrarily 
fixed and let k be such that 

(4.33) 

k < min {ki,k,2{\vo\, 16101), fi:3(||t;o||, |6'o|, Tq + r), K3(pi(r), \6o\,r)} 
■= fiii\\vo\\, \9o\), 

where Ki is given by (13.431) . /t2(-, ■) and K3(-, - , ■) are given in Proposi- 
tion 4-i. Tq is the time of entering the absorbing ball in H for {w", 6'"} 
and pi is given in fl4.37p below. Then we have 

(4.34) <K4(||^o||,|^o|), Vn>l, 

where K^^i^-, ■) is a continuous function defined on M^, increasing in 
both arguments. 
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Proof. In order to derive uniform bounds for ||f"||, for all n > 1, we 
apply Proposition 14.11 on successive intervals of time, with different 
initial values. On each interval considered, we obtain a bound 
that depends on the norm of the initial value (see (14.291) ). Using (the 
discrete uniform Gronwall) Lemma 14.41 we majorize the norm of each 
initial value by a constant pi (see (14.371) below) and recalling the fact 
that is an increasing function of its arguments, we obtain a bound 
independent of the initial value considered. 

We start by applying Proposition 14. li on the interval [0, To + r], where 
To := N^k, with A^^q being given in Corollary 13.21 We obtain: 



(4.35) 

and 
(4.36) 



\eo\,To + r),yn = 0,--- ,No + N, 
Nr:= [r/k\, 



V 



l + 2ciK^k 
Vn = I,-- - ,No + Nr. 
Rewriting (I4.36P in the form 

< ^n-l(l + kr]n-l) + kC, 



5z/ 



with in = \\v^W\ Vn = 2ciKl {\\v^f + IkM^) , (n = ^Mf, wc apply 
Lemma 14.41 with ni = A^o + 1, n2 = Nj. — 2, and n=K = A^^o + Nr. 
For n' = No + 1, No + 2, we compute, using (I3.3ip and (I4.33P and 
recalling that, by Corollary 13. 2[ the bounds Ki and Mi on and 
1 6*" I, respectively, can be replaced by po for n> No: 



n'+n2 



n'+n2 



rin = 2ciplk ^ 



,n||2 



kpl) 



2 2 

< -ClPo 



9 ^9 

Po + -Po^ 



n'+n2 



n'+n2 



14 2 . 14 2 



n=n' 

n'+n2 n'+n2 -, / , \ 

^E«" = *E ll""ll'si(pS + ipjir 

n— n' n— n' ^ 
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Then Lemma [4.41 implies: 

lu \r V J ou 

(4.37) / 2 4 A ^ 3 

exp < -CiPo I + 

=:pi(r)^ 

Using the above bound and recalhng assumption (14.331) on k and the 
fact that Ks{-,-,-) is a decreasing function of its arguments, we notice 
that k satisfies the constraint (14.141) in Proposition 14. 1[ We can, there- 
fore, apply Proposition 14.11 on the interval [Tq + r, Tq + 2r] with the 
initial data {v'^o+Nr ^^No+Nrj^ obtain that relation (iii) holds for all 
n = No + Nr + !,■■■ ,No + 2Nr, and 



\\v-\\<K,{\\v''^>+''^l\eo\,r)<KMr),\eolr) 
^ ' ' yn = No + Nr + lr-- ,No + 2Nr 



Using again Lemma 14.41 with rii = Nq + Nr + l,n2 = — 2 and 
n=K = iVo + 2Nj., we obtain 

(4.39) llt^^^^+'^' ll' <Pi(r)2. 

Iterating the above procedure, we find that < pi{r), for 

alH = 1, 2, ■ ■ ■ , and 

(4.40) < K,{p,{r), |^o|,r), Vn > iVo + N^. 

Finally, recalling (I4.35p . which gives a bound for < n < A^o + ^r, 
we obtain 

<max{i^3(||^°||,|^^o|,To + r),K3(pi(r),|^o|,r)} 

(4.41) 

=:^4(||t^o||,|^o|), Vn> 1. 
This completes the proof of the theorem. □ 
Corollary 4.1. Under the assumptions of Theorem \4.1\ we also have 

^ ll^n _ <Kl + tl^KlKlk{m -i + l) 

(4.42) n=i 

2 

+ -Mlk{m-i + l), \/i = 1,- ■ ■ ,m. 

Proof. Taking the sum of (14.111) with n from i to m and using (I4.34p 
gives conclusion (I4.42p of the corollary right away. □ 
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4.2. iJ^-Uniform Boundedness of 9^. We are now going to prove 
the if^-uniform boundedness of 6'", for all n > 0. In order to do so, we 
will first use the discrete Gronwall lemma to derive an upper bound on 
II ^"11, n < N, for some > 0, and then we will use another version of 
the discrete uniform Gronwall lemma (see Lemma 1^6] below) to obtain 
an upper bound on ||^"||, n > N. 

Lemma 4.5. Let {vq, 6q} G V and {f 6*"} be the solution of the nu- 
merical scheme (12.461) - fl2.47p . ^4/50, let T > be arbitrarily fixed and 
k be such that 

(4.43) k < min {k^HIvoI |^o|), f^^ihol |^o|)}, 



where K4(-, ■) is given in Theorem \4-l\ and 

1 27 c'^ 

(4.44) KM\v4,m = ™« = ^. 

Then we have 

(4.45) ||r f < A|^0||2 ^ J_ \ ^ vn = 1, ■ ■ ■ , iV := [T/k\. 
Proof. Replacing 6 by 2kA20"' in f l2.47p . we obtain 

(4.46) " " " " " " 21 , , 2 ; 

-2A;(t;^,A2r) + 2/tA;|A2^"|2 = 0. 

Using property f l2.38p of the trilinear form 62 and recalling fl3.3Up and 
fl4.34p . we have the following bound of the nonlinear term, 

2kb2{v'',9'',A2e'') < 2cfeA;|t;"|i/lt;"f/l^"||i/V2^"l^^^ 

(4.47) K 

< —i 

2 

Using the Cauchy-Schwarz inequality and recalling f l3.30p . we have the 
following bound 

-2A;(f 2", ^2^) < 2A;|fJP2r| 

- 2 ' 'ft ^ 
Relations (H^ -( H:i8|l imply 

||rf - \\9"-^f + ||r - ^-^2 ^ ^kiA^e^^i^ 

(4.49) 2 

< C2K!Klkm' + -K',k, 

K 

from which we obtain 

(4.50) lirf < -ll^^^i^ + —K?k, 

a K,a 



< -A;|A2r|2 + C2K^Klk\\9''f. 
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where 

(4.51) a = 1 - C2KlKlk. 
Using recursively (14.501) . we find 

(4.52) ||r f < (1 - C2KlKlk)-'' (pY + ^ 
Since 

l-x>4-'', 0<x<^, 

and, by hypothesis, C2KfKlk < 1/2, conclusion (14.451) follows imme- 
diately. This completes the proof of Lemma 14. 5[ □ 

In order to derive an upper bound on ||^"||, n > N, we will need 
the following version of the discrete uniform Gronwall lemma, slightly 
different from Lemma [4. 4t 

Lemma 4.6. We are given k > 0, positive integers no, rii and positive 
sequences ?7„, (n such that 

(4.53) kr]n < ^, forn > Uq, 

(4.54) (1 - kr]n)^n < U-i + kCn, for n > Uq. 
Assume also that 

fco+ni fco+ni 

k'^ r]n< ai(no,ni), ^ Cn < 02(^0,^1), 

(4.55) 

n=ko 

for any ko > Uq. We then have, 

(4.56) < (^^^^ + -2{no, nO) e^^^--), 

for any n > Uq + ni. 

Proof. Let and n4 be such that uq < n2 < < n2 + rii. Using 

recursively (14.541) . we derive 

(4.57) 
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Using the fact that 1 — x > e'"^^, Vx G (O, i), and recalhng assumptions 
(I4.55p ^ and M.SSp ^- we obtain 

Multiplying this inequality by k, summing ns from n2 + l to n2 + ni and 
using assumption M.SSP g gives the conclusion (14.561) of the lemma. □ 

We are now able to derive an upper bound on ||6'"||, n > N. More 
precisely, we have the following: 

Lemma 4.7. Let {vq, 6q} G V and {f 6'"} be the solution of the nu- 
merical scheme (12.461) - (12. 47p . Also, let T > be arbitrarily fixed and 
k be such that 

(4.58) k < mm\Ki{\\vo\\,\9o\),K5{\\vo\\,\eo\)J 



where /t4(-, ■) is given in Theorem \4.1\ and is given in Lemma 

\4-5\ Then there exists M2 = 71^2(1117011)1^01)^); given in (14.611) below, 
such that 

(4.59) liril <M2(||T;o||,|^^o|)T),Vn>iV:= [T/k\. 
Proof. We apply Lemma 14.61 to ( 14.49^ , which we rewrite as 

(4.60) 2 

We set a = lirf , 7]n = C2KIKI Cn = iKl no = 1, ni = iV - 1 and 
for fco > 1 we compute: 

kY,Vn = kY, C2KIKI < C2KIKIT, 
n=k(} n=ko 

fco+ni fco+ni n n 

kY.Cn = kY,'-Kl< '-KIT, 



K 

n=ko n=ko 



ko+ni ko+ni ^ j^2 \ 

^ E ^" = ^ E ll^"f <MMi^ + ^r (by(l332D). 



n=ko n=ko 

Then Lemma [4.61 implies 



||^n||2 + ^ + kit) e^^^^?^!^ 

(4.61) " " - V ^ / 

:=Mi{\\vol\Oo\,T),W>N. 

Thus, the lemma is proved. □ 
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Combining Lemma 14.51 and Lemma 14.71 we obtain that 6'" are uni- 
formly bounded in V, for all n > 0. More precisely, we have 

Theorem 4.2. Let {fo,6'o} ^ V <^f^d {"i^", 6*"} he the solution of the 
numerical scheme fl2.46p - fl2.47l) . Also, let T > be arbitrarily fixed 
and k be such that 

(4.62) k < min L^{\\vo\\, |^o|), ^^sdl^oll, l^^ol), ^1 =: /^edl^oll, l^ol), 



where /t4(-, ■) is given in Theorem \4.1\ and n^i^-,-) is given in Lemma 
\4-5\ Then there exists M3 = M3(||t>o||, ||^^o||); such that 

(4.63) ||^"|| <M3(||^;o||,||^^o||), Vr2>0. 
Proof. Taking 

Msdl^oll, ll^oll) = max |4^^^?^4^ (^||^of + ^) , M^iWvol |0o|,T)| , 

Lemmas 14.51 and 14.71 give conclusion f l4.63p of the theorem. □ 
Corollary 4.2. Under the assumptions of Theorem \4.S^ we also have 

m 

J^Wd""- ^""1^ <M^ + C2KfKjM^k{m - n + 1) 

(4.64) n=i 

2 

H — K^k{m — n + l), Vz = l,---,m. 

Proof. Taking the sum of (14.490 with n from i to m and using f l4.63p 
gives conclusion (14.641) of the corollary right away. □ 

Theorem 14.11 and Theorem 14.21 can be combined to obtain the fol- 
lowing 

Theorem 4.3. Let {^0,6*0} ^ V <^''^d {t>", 6'"} be the solution of the 
numerical scheme (12.460 - (12.470 . Then there exists a decreasing positive 
function k,7{-) and an increasing positive function K^l-), such that if 

(4.65) k<Ky{\\{vo,eo}\\), 
then 

(4.66) \\{v\e^}\\<K,{\\{vo,eo}\\),W>0. 
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5. Convergence of Attractors 



In this section we address the issue of the convergence of the attrac- 
tors generated by the discrete system fl2.45l) - fl2.47l) to the attractor 
generated by the continuous system fl2.lip - fl2.18p . Whereas for the 
continuous system fl2.lip - fl2.18p one can prove both the existence and 
uniqueness of the solution (see, e.g., [15])-and, therefore, define a global 
attractor-, for the discrete system fl2.45p - fl2.47p one can prove (using 
Theorem 14. 3 p the uniqueness of the solution provided that k < fi:(||Mo||), 
for some /€(||uo||) > 0. Since the time restriction depends on the initial 
data, one cannot define a single-valued attractor in the classical sense, 
and this is why we need to use the attractor theory for the so-called 
multi-valued mappings. Multi-valued dynamical systems have been in- 
vestigated by many authors (see, e.g., [1], [2], [1], [11], [12], [I3]), but in 
this article we use the tools developed in [5] to study the convergence of 
the discrete (multi- valued) attractors to the continuous (single- valued) 
attractor. For convenience, we recall those results in Subsection 15. H 
and then we apply them to the thermohydraulics equations in Subsec- 
tion O 

5.1. Attractors for multi-valued mappings. Throughout this sub- 
section, we consider (i/, | ■ |) to be a Hilbert space and T to be either 
M+ = [0, oo) or N. 

Definition 5.1. A one-parameter family of set-valued maps S{t) : 
2^ 2^ is a multi-valued semigroup (m-semigroup) if it satis- 
fies the following properties: 

(5.1) S{0) = I2H (identity m 2^); 

(5.2) sit + s) = S{t)S{s), for all t,seT. 

Moreover, the m-semigroup is said to be closed if S{t) is a closed 
map for every t G T, meaning that if Xn ^ x in H and yn G S{t)xn 
is such that yn ^ y in H , then y G S{t)x. (To simplify the notation, 
hereafter we have written S(t)x in place of S{t){x} .) 

Definition 5.2. The positive orbit of B, starting attET, is the set 



[jS{T)B, 



T>t 



where 



S{t)B 



U ^w^- 
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Definition 5.3. For any B E 2^ , the set 

is called the u-limit set of B. 

Definition 5.4. A nonempty set B E 2^ is invariant for S(t) if 

S{t)B = B, Vt G T. 

Definition 5.5. A set Bq G 2^ is an absorbing set for the m- 

semigroup S(t) if for every bounded set B E 2^ there exists ts E T 
such that 

S{t)B C Bo, Vt > te. 

Definition 5.6. A nonempty set C E 2^ is attracting if for every 
bounded set B we have 

lim dist{S{t)B,C) = 0, 

t—^co 

where dist(-, ■) is the Hausdorff semidistance, defined as 

(5.1) dist(i3,C) =supinf |6-c|,Vi3,C C H. 

beB 

Definition 5.7. A nonempty compact set A E 2^ is said to be the 
global attractor of S{t) if A is an invariant attracting set. 

Remark 5.1. The global attractor, if it exists, is necessarily unique. 
Moreover, it enjoys the following maximality and minimality properties: 

(i) if A is a bounded invariant set, then A D A; 

(ii) if A is a closed attracting set, then A E A. 

Definition 5.8. Given a bounded set B E 2^ , the Kuratowski mea- 
sure of noncompactness a{B) of B is defined as 

a{B) = inf |(5 : B has a finite cover by balls of X of diameter less than 5}. 
We note that a{B) = if and only if B is compact. 

The following theorem, whose proof can be found in [5] , gives condi- 
tions under which a global attractor exists. 

Theorem 5.1. Suppose that the closed m-semigroup S{t) possesses a 
bounded absorbing set Bq E 2^ and 

(5.2) lim a{S{t)Bo) = 0. 

t—^oo 

Then co{Bo) is the global attractor of S{t). 
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For the purpose of this article, we need to introduce the notion of 
discrete m-semigroups. More precisely, we have the following: 

Definition 5.9. Given a set-valued map S : 2^ ^ 2^ , we define a 
discrete m-semigroup by 

S{n) = Vn G N, 

and we will denote it by {S}nen (instead of {S'^}nef'i)- 

Remark 5.2. Given two nonempty sets B,C E 2^ , we write 

B-C = {b-c: beB,ceC} and |i3|=sup|6|. 

beB 

In order to prove the convergence of the attractors generated by 
the discrete system f l2.45p - fl2.47p to the attractor generated by the 
continuous system fl2.lip - fl2.18p we will use the following result, whose 
proof can be found in [5]; see also [21], [T9] . 

Theorem 5.2. Let S(t) be a closed m-semigroup, possessing the global 
attractor A, and for kq > 0, let {Sk, < k < KojneN be a family 
of discrete closed m-semigroups, with global attractor Ak- Assume the 
following: 

(HI) [Uniform boundedness] : there exists Ki G (0, Kq] such that the 
set 

}C= \J Ak 

A:e(0,/ti] 

is bounded in H ; 

(H2) [Finite time uniform convergence]: there exists tQ > such that 
for any T* > to, 

lim sup 15*^0; — S{nk)x\ = 0. 

''^^ x&Ak,nke[to,T*] 

Then 

limdist(A,^) = 0, 

where dist denotes the Hausdorff semidistance defined in (15. ip . 

5.2. Application: The thermohydraulics equations. The system 
fl2.lip - fl2.18p possesses a unique solution and thus generates a continu- 
ous single- valued dynamical system S(t) : H ^ H, with global attrac- 
tor A, bounded in V (see, e.g., [I5]). Using Theorem 14. 31 one can prove 
that the discrete system fl2.45p - fl2.47p has a unique solution provided 
that k < k{\\uo\\), for some /tdl-uoll) > 0. The dependence of the time 
step k on the initial data prevents us from defining a single-valued at- 
tractor in the classical sense, but this difficulty can be overcome by the 
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theory of the muhi- valued attractors. More precisely, in this article we 
will prove that there exists kq > such that if < k < kq, the system 
f l2.45p -( l2.47p generates a closed discrete m-semigroup {S'fc}„gN, with 
global attractors Ak, that will converge to A in the sense of Theorem 

[El 

In order to do that, we define, for > 0, the multi- valued map 
Sk '■ 2^ 2^ as follows: for every u = {v, 6} G H, 

SkU = {u = {v^O} & V : u solves fl5.3p - fl5.4p below with time-step k} : 
(5.3) 

(f , v') + vk{{v, v')) + kbi{v, V, v) — k{e29, v') = {v, v'), Wv' G Vi, 
(5.4) {6, 6') + Kk{{e, 6')) + khiv, e, 6') - k{v2, 6') = {6, 9'), W G V^. 

We then have the following: 

Theorem 5.3. The multi-valued map St associated with the implicit 
Euler scheme (12.451) - fl2.47p generates a closed discrete m-semigroup 

Proof. Since conditions (S.l) and (S.2) are satisfied by definition, we 
just need to prove that for each n G N, S'^ is a closed multi- valued 
map. For that, we let n G N be arbitrarily fixed and, as j — )■ oo, we let 
M° ^ M° in H, where = = {v^,9'^}. Also let u] G 

be such that m" — m" in H, where -u" = {f",6'"},M" = {t>",6'"}. We 
need to show that m" G SJ^u^. 

Indeed, since G SJ^u'^, there exists a sequence (u°, m], . . . , u^~^, u"^), 
with G SkU^~^, such that 

(5.5) 

iv],v') + ukiiv],v')) + khiv],v],v') - kie2e],v') = iv]-\v'), W G V,, 
(5.6) 

{9], 9') + Kk{{9], 9')) + kh2{v], 9], 9') - kiiv])2, 9') = (9]-', 9'), W G V2. 

The sequence m° being convergent in H, it is also bounded in H and 
thus there exists M > such that 

(5.7) sup|u°|2<M. 

j 

Then Lemmas 13.21 and 13.31 imply that for every i = 1, . . . , n, the se- 
quences Vj and 9j are bounded in Vi and V2, respectively. We therefore 
have that there exist subsequences still denoted Vj and 6**, such that as 
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j oo: 

(5.8) — )■ v\ strongly in Hi and weakly in Vi, 

(5.9) 6j — 6**, strongly in H2 and weakly in V2. 
Now, passing to the limit in fl5.5l) - fl5.6p . we obtain 
(5.10) 

{v\ v') + uk{{v\ v')) + kbi{v\ v\ v') - k{e29\ v') = {v'-\ v'), W e V^, 
(5.11) 

{e\ e') + Kk{{e\ e')) + kh2{v\ e\ e') - k{{v% e') = {e'-\ e'), w e V2. 

We therefore obtain that G SkU^~^, for each i = 1, . . . , n, and hence, 
G SkU^"^ C S]^vP . This completes the proof of the theorem. □ 

In order to prove the existence of the discrete global attractors, we 
first prove the existence of absorbing sets. More precisely, we have the 
following: 

Proposition 5.1. There exists Kg > 0, independent of {vo,dQ},n,k, 
such that if k & (0, Kg] the following holds: there exists a constant 
Ri > such that for every R > and \{vo,9o}\ < R, there exists 
Ni = Ni{R, k)>0 such that 

(5.12) \\Sj:{vo,eo}\\<Ru Vn>iVi. 
Hence, the set 

Bi = {{v,9}EV: \\{v,9}\\<Ri} 

is a V -hounded absorbing set for {Sk}n&h for k G (0, Kg]. 

Proof. Let ki be as in Corollary 3.2 and let k < min{l,Ki}. Also, 
let -R > and |{t;o,6'o}| < R. Then, by Corollary 3.2, there exists 
No = No{R, k)>0 such that 

(5.13) IK,r}|<po, yn>No. 



Let m := A^Vi + 



Then equations fl3.3ip and fl3.32p imply 



m 

(5.14) uk J2 ll^'T < Po + -Po(^-^o)A;, 



V 

j=No+l 



m ^ 

(5.15) nk \\e^f < pI + -pl{m- No)k. 



j=No+l 
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Adding the above relations we obtain 
(5.16) 

J2 Hv^\^ + K\\e^f)] < pI (2 + ^{m - No)k + ^{m - No)kY 

Assuming that for every j G {Ai'o + 1, ■ ■ ■ , m} 
{v\\v^f + Kll^^'f ) > ^° (2 + -(m - Aro)fc + i(m - Aro)fc 



/c(m — iVo) \ z/ K 

we obtain 
(5.17)^ 

k{ ^ (z/||t;i2 + /t||^i2) 1 >p^2 + i(m-iVo)A; + -(m-iVo)A;y 

which contradicts (15.161) . Hence there exists / G {Nq + 1, ■ ■ ■ , m} such 

that 

(5.18) 

{y\\v'f + K\\&f) < , , ^° ^ , (^2 + -{m - No)k + -(m - No)k 
k[m — i\q) \ u k 

We, therefore, have 

(5.19) \\{v\9^}r<2pl(2 + ^ + ^) fl + i) =:Rl 

Applying Theorem 14.31 with initial data we obtain that there 

exists K^{\{v\d^}\) and K^{\{v\Q^}\) such that if < k^{\{v\B^}\), 
then 

(5.20) \\{v\ff^}\\<K,{\{v\&}\),^n>l. 

Recalling (I5.19P and the fact that K7(-) and K^{^ are, respectively, 
decreasing and increasing functions of their arguments, (I5.20p yields 

(5.21) ||K,^^"}II < ^5(^*) =: i?i,Vn > A^i = N^{R,k) := N^^ 
provided that k < Kg, where 

(5.22) Kg = min{l, Ki, K7(_R*)}. 

This completes the proof of Proposition 15.11 □ 

We are now in a position to prove the existence of the discrete global 
attractors. More precisely, we have the following: 
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Proposition 5.2. For every k G (0, ks], there exists the global attractor 
Ak of the m-semigroup {Sk}neN- 

Proof. Let Bq = Bh{0, po) be the bounded absorbing set given in Corol- 
lary [3]2l Then Proposition 15.11 implies that SJ^Bq is bounded in V, for 
all n > Ni{pQ,k). Since V is compactly embedded in H, we obtain 
that SJ^Bq is relatively compact in H and, thus, a{SJ^Bo) = 0, for all 
n > Ni{pQ,k). Condition f l5.2p of Theorem 15.11 is therefore satisfied 
and then the existence of the discrete global attractor Ak follows right 
away. □ 

Remark 5.3. Since the global attractor Ak is the smallest closed at- 
tracting set of H, Proposition \5. 1\ implies 

(5.23) Aci3i,VA;G (0,fi:8], 
and thus 

(5.24) U AkdBi. 

fce(o,K8] 

Let us recall that our goal is to prove, using Theorem 15. 2[ that the 
discrete global attractors Ak converge to the continuous global attrac- 
tor A. Thanks to fl5.24p . condition (HI) of Theorem 15.21 holds true. 
There remains to prove the finite time uniform convergence required 
by (H2). In order to do that, we define, for any /c > and for any 
function ip, the following: 

(5.25) 7/;fc(t)=^", te[{n-l)k,nk), 

(5.26) ^,(t)=r + ^^{r-r-'), te[{n-l)k,nk). 

With the above notations, equations (12.461) and fl2.47p can be rewrit- 
ten as follows; for t G [(n — l)k, nk): 

(5.27) 
f dvkjt) 
V dt - 
(5.28) 

'd Okit) 
dt 

where 
(5.29) 



V ) +u{{vk{t),v)) + bi{vk{t),Vk{t),v) = {e29k{t),v) + Uk{t),v), Wv G V^, 

+ KiiOkit), 6)) + b2{vk{t), Okit), 6) - {vk{t))2, e) = {gk{t), e),\/ee V2, 

{fk{t),v) = v{{vk{t) - Vkit),v)) + bi{vk{t), Vk{t), v) 

-bi{vk{t),Vk{t),v) - {e2{9k{t)-ek{t)),v), 
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-b2{vk{t),ek{t),e)-{{dk{t)-vk{t))2,e). 

Lemma 5.1. Let T* > be arbitrarily fixed and let k < Kq, where 

(5.31) Ko = min{K8, KjiRi)}, 

with Kg being given in fl5.22p and n-j being given in Theorem \4-^ As- 
sume that {fo, Oq} E Ak and let {f " , be the solution of the numerical 
scheme (EIS]) -(123ZD ■ Then there exist K^{T*) and K71t*) such that 

(5.32) \\fk\\l2^o,T*;Vi) < kMT*), 
and 

(5.33) lkfc|li2(o,T*;V2') < kKj{T*). 

Proof. Let us first note that for any t E[{n — l)k,nk) we have 

(5.34) V^,(t) - ^,(t) = ^^(^- - r-'). 

k 

Also, since {^0,6*0} G Ak, we have that ||{fo,^o}|| ^ Ri (by fl5.23p ) 
and then Theorem 14.31 imphes that for k < kq, 

(5.35) \\{v\9^^}\\<K,{R,),\fn>0. 

Now let f G Vi be such that ||f || < 1, and let t G [(n — l)k,nk) be 
fixed. Using property fl2.3ip of the trilinear form hi, we have 

\bi{vk{t),Vk{t),v) - bi{vkit),Vkit),v)\ 

= \bi{vk{t) - Vk{t),Vk{t),v) + bi{vk(t),Vk(t) - Vk{t),v)\ 

<c.(||^^.(t)-^;.(t)||(||5.(t)|| + ||^;,(t)||)||t;|| 

<c||t;"-w"-l fbv and \\v\\ < 1). 

We also have 

(5.37) i^\{{vk{t) - Vk{t),v))\ < u\\v^ - v^-'\\, 

(5.38) \{e2{9k{t)-9k{t)),v)\ < 11^-^-1. 
Relations fO^ - dE^ imply 

(5.39) imnvi < edit;" - v--'\\ + iir - r-i), 

and thus, setting A^* = [T*/k\ and recalling that 0,6*0)11 < -Ri , we 
obtain 

-T* N*+l „„fc 



(5.40) 



\\fk\\h(o,T-vi)= ll/fcWlly/rfi= E / \\Mt)\\l,dt 

Jo J{n-l)k 



<kK,{T*) (by (15311), (!132D,(!M1), 
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which proves f l5.32p . 

Now let 6 e V2he such that ||6'|| < 1, and let t e [{n - l)k,nk) be 
fixed. Using property fl2.36p of the trilinear form b2, we have 

(5.41) 

Mvk{t)Jk{t),9)-h{vk{t),9k{t),9)\ 

= I b2 (vk (t) - Vk {t)Mi).d) + h2 {vk {t)Mt)-Ok{t),e)\ 

< ctiiMt) - Vkimwokim + Hmmt) - e^mmi 

< c(||t;" - v^'-^W + ||r - r^i) (by dESl, (15351) and ||^|| < 1). 
We also have 

(5.42) f^imit) - 9kit),e))\ < k\\9'' - 9''-'\\, 

(5.43) \iidkit)-Vkit))2,9)\ < \\v^-v^-\ 
Relations (l53I])-(IE43]) imply 

(5.44) \\9kmvi < ci\K - + r-i), 

and thus setting N* = [T*/k\ and recalling that ||{t'0)^o}|| < -Ri , we 
obtain 

,^ lkfc|li2(o,T*;y') = / \\9k{t)\\v'dt = / \\gk{t)\\l,dt 

(5.45) ' ' ^' Jo ' ^ J{n-l)k ' 

<kKj{T*) (by (l521,ga2D,(|4S4D), 

which proves fl5.33p and the proof of the lemma is complete. □ 

We are now able to prove that condition (H2) of Theorem 15.21 is 
satisfied. More precisely, we have the following 

Proposition 5.3 (Finite time uniform convergence). For any T* > 

we have 

(5.46) lim sup \Sl{vo,9o} - S{nk){vo,9o}\ = Q. 

{yo,eo}e^fc,nfce[0,T*] 

Proof. Let 

(5.47) Ut) = <t) - Vk{t), Vk{t) = 9it) - 9k{t). 

Subtracting and ([52SD from (imj) and flTT^ written in their 

week form, respectively, we obtain 



(5.48) V ^ 

+ hi{vk{t),ik{t),v') = {e2r]k{t),v') - {fk{t),v'), W G Vi, 
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(5 49) + ^'^^ + ^'^ 

Replacing v' by C,kit) in (15.481) . we find 

Using property fl2.3ip of the form bi, we bound the nonhnear term as 

bi{ut),v{t),ut))<cbMmMm\\vm 

o u 
Using the Cauchy-Schwarz inequahty, we also have 

\ie2Vkit),um < Mt)\Mt)\ 

(5.52) ^ \VkmMt)\\ 

<^ii6(t)ir+-k(t)r, 

mt),um<\m)\\vi\Mt)\\ 
^^■^^^ <"7:\Mt)r + -\\hmh. 

6 u 1 

Relations imply 

;^ie.(t)r+Hia(t)ir<-ii^(t)fia(t)p 

(5.54) dt V 

Now replacing d' by ri^^if) in (15.491) . we find 

-((a(t))2,r/fe(t)) = -((7fc(t),r7fc(t)). 
Using property fl2.36p of the form 62, we bound the nonlinear term as 
(5.56) 

|fe2(e.(t),^(t),r/.(t))|<c,|e,(t)r/2||^,(t)ir/^||^(t)|||r^,(t)|V2||^,(t)||V2 

^l\\mf\m?^\\mfw{t)?- 
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Using the Cauchy-Schwarz inequality, we also have the following bounds: 

mk{t))2,vkm<Mmvk{t)\ 

O K 



\{9k{t),Vk{t))\ < \\9k{t)\\v^\\Vk 

^'■''^ <-jvkit)r+-\Mt)rn- 

Relations fl535ll -f l538l) imply 

j^\vk{t)\' + K\\r^k{t)r ^Ukitw + ^mt)rMt)\' 

(5-59) + -\\0{t)r\Vk{t)\' + -Mt)\' 

K K 



C , 



2 



(5.60) 



\\9kK'^J\\v^ 

Adding equations (15.541) and (I5.59p . we obtain 

|(ie.(t)P + k(t)P) + ^Hle(t)f + «:h(t)|P 

As shown in pLSj, the solution {f , ^} of the continuous problem is uni- 
formly bounded in V for all t > 0. More precisely, we have 

(5.61) sup sup \\S{t){v^M\\<(^- 
Thus, inequality (I5.60p becomes 

^(i6(t)r + ir/.(t)r) + |Hie(t)ir + «:h(t)ir 

(5.62) 6 

< ci\m\' + Mm + ^Ukimh + ^iMmvi- 

By Gronwall's lemma and using the fact that .^^(0) = ri{0) = 0, we 
obtain 

(5.63) Mt)\' + \Vk{t)\'<ce^^\\\M 
and recalling (15.321) and (15.331) . we find 

(5.64) \Ut)\' + \Vk{t)\'' <ck, 
for some constant c = c(T*) > 0. 
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We therefore have, 
(5.65) 

hm sup \SJ^{vo,9o} - S{nk){vo,9o}\ 

= hm sup sup \{v'',6"-} - {v{nk),6{nk)}\ 

= hm sup sup \{vk{nk),9k{nk)} — {v{nk),9{nk)}\ 

{vo,eo}(^-Ak,nke[o,T'-] {v^,e"}esj!{vo,eo} 

= hm sup sup \{^k{nk),7]k{nk)}\ = 0, 

which concludes the proof of the lemma. □ 

We have, therefore, proved that conditions (HI) and (H2) of Theorem 
15. 21 are both satisfied and thus, the long-term behavior of the semigroup 
S(t) generated by the continuous thermohydraulics equations (12.111) - 
(I2.12P is approximated by that of the m-semigroups generated by the 
discrete system (I2.45p - (l2.47p . More precisely, we have the following: 

Theorem 5.4. The family of attractors {Ak}ke{o,Ko] converges, as k ^ 
0, to A, in the following sense: 

limdist(A,^) = 0, 

where dist denotes the Hausdorff semidistance in H , namely 

dist(^fc,^) = sup inf \xk — x\. 
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